home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / DPPDI.z / DPPDI
Encoding:
Text File  |  1998-10-30  |  1.8 KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPPPPPDDDDIIII((((3333FFFF))))                                                            DDDDPPPPPPPPDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPPDI   - DPPDI computes the determinant and inverse of a double
  10.      precision symmetric positive definite matrix using the factors computed
  11.      by DPPCO or DPPFA .
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE DPPDI(AP,N,DET,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAAPPPP DOUBLE PRECISION (N*(N+1)/2)
  21.         the output from DPPCO or DPPFA.
  22.  
  23.      NNNN INTEGER
  24.         the order of the matrix  A .
  25.  
  26.      JJJJOOOOBBBB INTEGER
  27.         = 11   both determinant and inverse.
  28.         = 01   inverse only.
  29.         = 10   determinant only.  On Return
  30.  
  31.      AAAAPPPP the upper triangular half of the inverse .
  32.         The strict lower triangle is unaltered.
  33.  
  34.      DDDDEEEETTTT DOUBLE PRECISION(2)
  35.         determinant of original matrix if requested.
  36.         Otherwise not referenced.
  37.         DETERMINANT = DET(1) * 10.0**DET(2)
  38.         with  1.0 .LE. DET(1) .LT. 10.0
  39.         or  DET(1) .EQ. 0.0 .  Error Condition
  40.  
  41.      AAAA division by zero will occur if the input factor contains a zero on the
  42.      diagonal and the inverse is requested.  It will not occur if the
  43.      subroutines are called correctly and if DPOCO or DPOFA has set INFO .EQ.
  44.      0 .  LINPACK.  This version dated 08/14/78 .  Cleve Moler, University of
  45.      New Mexico, Argonne National Lab.  Subroutines and Functions BLAS
  46.      DAXPY,DSCAL Fortran MOD
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.